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Abstract 

— j The effective action for baryons in the external electromagnetic 

fields is obtained on the basis of the QCD string theory. The area law 

C^ I of the Wilson loop integral, the approximation of the Nambu-Goto 

,_±- ' straight-line strings, and the asymmetric quark-diquark structure of 

nucleons are used to simplify the problem. The spin-orbit and spin- 

*Q ', spin interactions of quarks are treated as a perturbation. With the 

help of the virial theorem we estimate the mean radii of nucleons in 

ON ■ terms of the string tension and the Airy function zeros. On the basis 



ON 



X 



of the perturbation theory in small external electromagnetic fields we 



i-C I derive the electromagnetic polarizabilities of nucleons. The electric 

11 ■ and diamagnetic polarizabilities of a proton are a p = 10 x 10 -4 fm 3 , 

Ph. (5p ta = — 8 x 1CT 4 fm 3 and for a neutron we find a n = 4.2 x 10~ 4 fm 3 , 

c| , (3% m = —5.4 x 10 -4 fm 3 . Reasonable values of the magnetic polar- 

izabilities J3 P = (5 ± 3) x 10" 4 fm 3 , j3 n = (7.6 ± 3) x 10" 4 fm 3 are 
estimated using the A contribution to the paramagnetic polar izability 
of the nucleons. 



1 Introduction 

In recent papers [1,2] we have calculated some electromagnetic characteristics 
of mesons in the framework of the QCD string theory [3,4]. Starting with 
QCD and excepting the area low for the Wilson loop integral, the effective 
action for mesons in external electromagnetic fields was derived. This ap- 
proach takes into account the chiral symmetry breaking and the confinement 
of quarks [3,4]. This work is the continuation of [5] on the case of baryons, 
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i.e. the three-quark system. We make some assumptions such as (i) the 
area law behavior of the Wilson loop, (ii) the straight-line approximation 
of the strings, (iii) the handling of spin degrees of freedom as a perturba- 
tion, (iiii) the quark-diquark structure of baryons. The spin-orbit, spin-spin 
interactions and Coulomb like short-range contributions are also neglected, 
as a first approximation. The area-law of the Wilson loop integral (or the 
confinement of quarks) was confirmed by the Monte-Carlo simulations and 
by the theoretical considerations. One of the approaches which explains the 
confinement of quarks is the approach [6] introducing stochastic gluon vac- 
uum fields. The second one uses the dual Meissner effect to guarantee the 
confinement of colour [7,8]. In our consideration, the confinement of quarks 
is postulated by excepting the area-law of the Wilson loop. The use of the 
straight-line approximation for the Nambu-Goto string is useful to simplify 
the calculations. Neglecting short-range interactions is justified by consider- 
ing large distances. The quark-diquark structure of baryons was considered 
in the framework of the nonrelativistic [9] and relativistic [10] models. In this 
case the linear baryon Regge trajectories have the same slope as for mesons 
[3]. 

The paper is organized as follows. In Section 2 we describe the gen- 
eral background and derive the Green function of a baryon. The effective 
action for baryons in external electromagnetic fields based on the proper 
time method and Feynman path-integrals is found in Section 3. Section 4 
contains the calculation of average distances between quarks and electric po- 
larizabilities of nucleons. We derive diamagnetic polarizabilities of protons 
and neutrons using the perturbative expansion in the small magnetic fields 
in Section 5. Section 6 contains a conclusion. 

2 The Green function of three-quark system 

Here we derive the Green function of three-quark system using the Schwinger 
proper time method and the Feynman path-integrals. Our goal is to calcu- 
late some electromagnetic characteristics of nucleons. For this purpose we 
need the effective action for baryons in external electromagnetic fields. The 
method of the Green functions will be explored. Let us consider the Lorentz- 
covariant and gauge invariant combination of three-quark, colourless system 
(baryon) [4] 



X B (x, y, z, d) = e abc [$(Z„, x)q(x)] a [$(Z , y)q(y)] b [$(Z , z)q(z)} 



c ' 



where q(x) is a quark bispinor; a, b, c are colour indexes so that [&(Z , x )l{ x )] a 
$ aa ' (Z , x)q a '(x) and e a b c is the Levi-Civita symbol (ei 2 3 = 1). We imply that 
quark fields q(x), q(y) and q(z) possess the definite flavours which will be 
specified later at the considering nucleons. As usual there is a summation on 
repeating indexes. The gauge invariance is guaranteed here by introducing 
the parallel transporter [4]: 



$(Z„, x)=P exp { ig / Apdzy. , (2) 




where P is the ordering operator along the contour C\ of integration, g is 
the coupling constant, A^ = A a \ a \ A a are the gluonic fields; A a are the 
Gell-Mann matrices and Z is the string junction. Although we have an 
attractive feature - gauge invariance, the function (1) depends on the form 
of the contour C. As X B (x, y, z, Cj) is a gauge invariant object, it obeys the 
Gauss law on the spacelike surface E. Three quarks (with fields q(x), q(y) 
and q(z)) are situated in the four-points x, y, z, respectively and connected 
with the four-point Z which is arbitrary four-point. Later, the position of Z 
will be defined by requiring to have the minimal area for the world surface of 
three-quark system [4]. The path of integration in Eq. (2) is also arbitrary. 
Let the four-points x, y,z and x', y',z' be the initial and final positions of 
three quarks, respectively. Two particle quantum Green function is defined 
as [4]: 

G(xyz,x'y'z') = {X B (x,y,z,C i )X B (x',y',z',Cl)), (3) 

were X B (x', y', z', C'A corresponds to the final state of a baryon: 

X B (x' : y\ z\ C>) = e mnk [q(x'Mx' : Z' )] m WW, Z' )] n [q(z')$(z' : Z' )] k . 

(4) 
Here q = q +r y4, q + is the Hermitean conjugate quark field; 7^ are the Dirac 
matrices and the brackets (...) mean the path- integrating over gluonic and 
quark fields: 

(X B X B ) = J DqDqDA^ exp {iS QCD } X B X B , (5) 

with the QCD action Sqcd- We imply that the measure DA^ in the path 
integral (5) includes the well known weight for gluonic fields [11]. The 
Minkowski space is used here but it is not difficult to go into Euclidean 
space to have well defined path-integrals. 



Let us introduce the generating functional for Green's function to calcu- 
late the path-integral (5) with respect to quark fields: 

m „] = I Dm exp {iS QCD + ij dx ( i(l K(.) + *<*)*<*))} , (6) 

where we introduce the external colour anticommutative sources r] a , fj a . Then 
the Green function (3) can be written as 



eabcemnk$aa>(Z ,x)$w(Zo,y)®cc>(Z ,z) 



G(xyz,x'y'z') = J DA^ 

x$ m / m (a/, Z' Q )$ n i n (y', Z' )<S> k , k (z' ', Z' ) 

x 5 6 /5f] af (x)5f] b >(y)5f] cf (z)5r] ml (x')5r] nr (y')5r] kr (z')Z[fj,rj} 



• (7) 

r/=r/=0 

In Eq. (7) Zq and Z' are the initial and final positions of the string junc- 
tion, respectively. The total surface which consists of world motions of quarks 
and the path of the string junction must be minimal. This requirement de- 
fines the path of the string junction [4]. Now it is possible to integrate the 
path-integral in Eq. (7) over quark fields q, q as expression (6) is a Gaussian 
integral. We may represent the QCD action in the form 

Sqcd = S(A) - I dxq(x) (7 A + m) q(x), (8) 

where S(A) is an action for gluonic fields with the included ghost fields, 
D^ — 9^ — igA^; m is the quark mass matrix and we imply the summation 
on colour and flavour indexes. Inserting Eq. (8) into Eq. (6) and integrating 
with respect to quark fields, we arrive at the expression 

Z[fj, 77] = det(-7 M £> M - m) exp I iS(A) + i / dxdyfj(x)S(x, y)rj(y) > , (9) 

where the classical quark Green function S(x, y) is the solution of the equa- 
tion 

ij tJL D ll + m)S(x J y) = S(x-y). (10) 

Using Eq. (9) and calculating the variation derivatives in Eq. (7) we find 
the quantum Green function of a baryon: 



G(xyz, x'y'z') = / DA^ det(-7 M £> M - m) exp {iS(A)} 



^abc^mnk 



X ^>bm\yi X )^an\ X iy )^ck\ Z y Z ) ^am\ X i X J^bniV 'iV ')^ck\ Z y Z ) 

+ S cm {z, x )b bn (y, y )S ak (x, z ) — S cm (z, x )S an (x, y )b bk (y, z) (11) 

+\m( 3; ! X )^cn\ Z i V J^bkiy^ Z ) ~~ ^fernl?/' X )^cn\ Z i 2/ )^ak\ X i Z ) ) 

where we introduce the following notation for the covariant Green function 

S am( X i X ') = ®aa>(Zo,x)S alm >(x,x')® rn/rn (x', Z' Q ). (12) 

The functional determinant in Eq. (9) describes the contribution from the 
vacuum polarization and gives the addition quark loops. As a first approx- 
imation we neglect the contribution of the loops. The presence of different 
terms in Eq. (11) is connected with the permutations of quark fields because 
the quantum Green function being considered. As different terms in Eq. (11) 
have the same structure, we consider in detail only one term. Neglecting the 
functional determinant we find the approximate expression for the baryon 
Green function 

Gi(xyz,x'y'z') = -J DA^exp {iS(A)} e abc e mnk Sf m (x, x')S bn (y,y')Sf k (z, z'). 

(13) 
Expression (13) is the basic formula for deriving effective action for baryons. 

3 Effective action for baryons 

To calculate some electromagnetic characteristics of nucleons we consider 
baryons in external electromagnetic fields. Then covariant derivatives be- 
come -D M = <9 M — igAp — iQA e J, where A^ is the vector-potential of electro- 
magnetic fields, Q is the charge matrix of quarks, Q =diag(ei, e 2 , ..., e^), e.% 
are charges of quarks, Nf is the number of flavours. Using the proper time 
method and Feynman path-integrals we found in [2] the Green function of a 
quark which is the solution to Eq. (10): 



roo rz(s)=x / ^ 

S(x, x')=i ds Dz [m- -iu,zJt) 

JO Jz(0)=x' \ 2 



x Ps exp < i / dt 



-zlit) -m 2 + ez,(t)A^(z) + E^ (ei# + gF t 



$(x, x') 

(14) 



where F%, = d^A e J - d u A e ^ F^ = d^A v - d u A )l - ig [A^ A v \ are the strength 
tensors of electromagnetic and gluonic fields, respectively; E^ = — (i/4) [7^,7^] 
are the spin matrices; Ps is the ordering operator of the spin matrices E MJ/ ; 
Zuit) = dz^(t)/dt; Zuit) is the path of the quark with the boundary condi- 
tions 2^(0) = x'^, z^s) = x^ and $(x, x') is the path ordered product (see 
Eq. (2)). Inserting Eq. (14) into Eq. (12) we find 



O nrY . I X , Jb 



roo rz(s)=x / ^ 

-/ / ds I Dzlm- -l^z^t) ) P s 







z(0)=x' 



xexp \i 1 dt 



hl(t) - m 2 + ez»{t)A*{z) + £^(eP^ + gF, v ) \ \ (Q Cm (Z , Z' )) am , 



(15) 

where the contour C x in Eq. (15) consists of lines between Zq, x and Z' 0) x' 
and path z^it). Using the expression (15) for each quark, the baryon Green 
function (13) becomes 



m i-lp^\tj) 



dixyz, x'y'z') = % f°° J[ d Sj f JJ Dz ij) 
3 3 

xP s ]~Jexpi I' 3 dtjE 



'"Sa&itj) 



ex P yH jf dtj [\ {^\t 3 )f - m) + eizf{t 3 )A*(z®) + e-L^ 



x (W(C x C y C z )), 



(16) 



where j = 1, 2, 3; tj is the charge of the j'-th quark; the boundary conditions 

4 1} (o) = ^, zj»( 8l ) = *„, 4 2 )(o) = y ;, 4 2 )( S2 ) = y„, 4 3 )(o) = 4, 4 3 )( S3 ) = 

z M are used here and the Wilson loop is given by (see [4]) 

(W(C x C y C z )) = e abc e mnk (($ Cx (Z , Z' )) ($ Cy (Z (h Z / )) b ($ Cz (Z , Z' )) ck ). 

(17) 
We took into account that [4] 

exp [ig J* dffi^} W(C) = exp j/J dt^j-^ 



W(C), 



where a^ u {t) is the surface around the point z^it). 



Contours C x , C y , C z correspond to three quarks which have paths z^\ 
z< u\ z< u^ an d masses m 1; m 2 , rn 3 , respectively. Relationship (16) is the gener- 
alization of one [3,4] for the case of quarks placed in external electromagnetic 
fields which possess spins. 

We imply further that the average distance between quarks (r) is greater 
than the time fluctuations (in units c = h — 1) of the gluonic fields T g : 
(r) > T g . As T g ~ 0.2 -e- 0.3 fm [12,13] so the condition (r) > T g is valid not 
only for asymptotic baryon states of the Regge trajectories with large angular 
momenta of the baryon but also for lower baryon states. The asymptotic of 
the average Wilson loop integral obeys then the area law and is given by (in 
the Minkowski space): 

(W(C x C y C z )) = exp {-ia (S 1 + S 2 + S 3 )} , (18) 

where a is the string tension and Sj (j = 1, 2, 3) is the minimal surface 
bounded by the trajectories of the quark q^ and string junction Z . The path 
of the string junction is defined by the requirement that the sum S± + S 2 + S 3 
is minimal [3,4]. Following [3,4], new variables are introduced: 

t{T t 2 T t 3 T T 

T = = = , Vj = 7Ti ( 19 ) 

s l s 2 s 3 ^ s j 

where r means the proper time for every quark and fij (j = 1,2,3 ) is the 
dynamical mass of the j— th quark. Using Eqs. (18), (19) from Eq. (16) we 
arrive at 



G^xyz^'y'z') = -i- [°° U^f [l[Dz® (m, - W ^V)) x 

JO ■ fij J j 

xP s ]Jexpi — £^ / dr 
Y 1 2/ij Jo 

exp{i / (IrV 
Jo ; 



5auu[T) 



r (i) 

,2 




~ 2 N (iJV)) - ^ + e^(r)<(^)) + ^ V F* V 



-ia(S 1 + S 2 + S 3 )}. (20) 

The integral in the last exponential factor of Eq. (20) represents the 
effective action for three-quark system (baryon) taking into account the spins 
of quarks. So T is the time of the observation, r is the proper time of 



quarks; rrij and fij are the current and dynamical masses of j— th quark. It 
follows from Eq. (20) that there is integration over dynamical masses fij to 
get the Green function of the three-quark system. Pre-exponential factors 
in Eq. (20) allow us to calculate in principle the spin-spin and spin-orbit 
contributions to the effective action. As a first approximation (see [3,4]) we 
neglect the short-range spin corrections and consider therefore scalar quarks. 
The terms e^E^F^, which describe the interaction of the magnetic field with 
the spins of quarks will be omitted. With this assumption we arrive at the 
effective action for baryons in external electromagnetic fields 



B = i>? 



1 / ,.,, s \2 m 2 



■cr(Si + S 2 + S 3 



H^Hr)) --± + ej zj>\T)Ai(z<») 

Z/Ij 

(21) 
The case when electromagnetic fields are absent was considered in [3,4]. 
The terms which describe the interaction of quarks with electromagnetic 
fields are essential for us because we are going to calculate electromagnetic 
characteristics of nucleons. It is convenient to introduce new variables R^, 
£ M and i]^ instead of z^ in accordance with relationships [3,4]: 

M M \m (^ + /i 2 ) ) ^ U(^+/W " 

/ \ V2 / \ 1/2 

z m = r + ( &* ] e + ( m i „ 



z 




f) = R ^_mm±M\ & , (22) 



where i? M is the center of mass coordinate of a baryon; £ M and f]^ are relative 
coordinates of quarks, M = fi\ + \x 2 + 7/3 is the sum of dynamical masses 
of quarks. The arbitrary mass parameter \x in Eq. (22) defines the scale 
of relative coordinates £ M and rj^. After the substitution (22), the measure 
YljDz^ transforms into DRDr/Dt; in path-integral (20). 

Let us consider the uniform and constant external electromagnetic fields. 
Then the vector-potential of electromagnetic fields can be represented as 

Al l (z^) = \F^z\p. (23) 



Inserting Eqs. (22), (23) into Eq. (21) we find the effective action for the 
three quark system (see [5]) in the form 



B= f dr 

Jo 



1 



2 " 2 V" "V Y 2/i, 



^(Si + ^ + ^ + AB, (24) 



AB = -F^I dr 



+ 7 (^6 + 4^) + 5 (^ + %Rv) + & 






where we introduce parameters: 



7 



(25) 




A 



JL 

M 



Ml + /i 2 



(ei + e 2 ) /i 3 e 3 (/ii + /x 2 ) 



(26) 

Ml + M2 M3 

and e = ei + e 2 + e 3 is the charge of a baryon. As a particular case, when 
electromagnetic fields are absent (AB = 0) we arrive at the action derived 
in [3]. It follows from Eq. (24) that the center of mass coordinate R^ is 
separated from relative coordinates £^ and rj^ and \x plays the role of the 
mass of the £ M , 77^ excitations. Following [3,4], the straight line approximation 
for strings and the asymmetric quark-diquark structure of baryons will be 
assumed. The asymmetric configuration (see also [9,10]) means that two 
quarks q^ and q^ are near each other and quark q^ is farther from them. 
This case is preferable [9,10] and the slope of linear baryon Regge trajectories 

'ii + It + a) and 



is the same as for mesons [3,4]. Then |£|^H?7|(|£| : 

the coordinate 77^ can be ignored. We neglect therefore the surfaces Si, S 2 

and assume for S3 the following expression [3,4]: 



So 



/ dr U , 
Jo 



(27) 



where 

I H (/ii + /i 2 ) / /J/X3 / 98 n 

Equation (27) takes into account the confinement of quarks and gives the 
linear potential between quarks. Using the definition B — J drC, where L 
is the Lagrangian, from Eq. (25) we arrive at the effective Lagrangian for 
baryons 

£ = ^LrI + »g _ V Hi _ ab | £ | +£ eZ_ (2Q) 

Here we neglect the coordinate i]^ and introduce the notation: 
1 



2 "" 



eR^R u + A^^ + 7 (i?^ + £» R v)} ■ ( 3 °) 



Lagrangian £ eZ describes the electromagnetic interaction of the string. It 
follows from Eq. (3.22) that 



,(3) _ I (M) , M)\ ^ , 1 / /i (H2-li\ 



As the second term in Eq. (31) is small, the coordinate £^ is proportional 
to the "distance" between quark q^ and the center of mass of quarks q^ 
and q^ which form a diquark. At the large time T limit £4 = 0, R4 = ir 
[3,4] and Lagrangian (30) takes the form 

C e = e (RE) + 7 (£E) — -e mnk H k \eR m R n + A£ m £ n + 7(-R m £ n + £ m -R n )J , 

(32) 
where the electric field Ej, = iF k i and magnetic field H k = (l/2)ek mn F mn . 
To clarify the physical meaning of the terms in Eq. (32), let us consider the 
dipole moment of quarks. Using the definition of the electric dipole moment 
and Eq. (22) we have 

d = E e i z ° ) = eR + 7^ + ^- ( 33 ) 



So, first two terms in Eq. (32) ( neglecting the coordinate rj^ ) describe 
the interaction of the dipole moment of quarks with the electric field in 



10 



accordance with the expression for the potential energy: U = — (dE). The 
magnetic moment is given by 

— 1 V^ (J) • U) — 1 

iilk r\ rank 2_^, j m n Q^-mnk 



GJ^rn^n "T -^SmSra T7 I ^"mK,n "T sm-'Wi 



+ P?? m ?7n + 8 (RmVn + VmRn) + SJ- 3 , (£mVn + Vmt 



(34) 



It follows from Eqs. (25), (32) that there is an interaction of the magnetic 
field with the magnetic moment in such a way that the interaction energy is 
U = — (mH). So Lagrangian (32) describes the interaction of electric and 
magnetic moments of baryons with electric and magnetic fields, respectively. 

4 Mean relative coordinate and electric po- 
larizabilities of nucleons 

Now we estimate the mean size of baryons. At the large time T limit £ = 
(£4 = i£o), R& = iT [3,4] and therefore only three dimensional quantities R^ 
and £k are important. From Eq. (32) we find three momenta corresponding 
to the center of mass coordinate Rk and relative coordinate £& : 

n*; = - = MR k — -e knm H m (eR n + 27^ n ) , 
oR k 2 

dC eff • 1 

TTfc = ^ = fj£k - -EknmHm (A£ n + 2"/R n ) , (35) 

d£k 2 

Here we take into account that Lagrangian is defined within the accuracy 
of the total derivative on time. Then the effective Hamiltonian Ti e ff = 
H-kRk + T^k^k — £eff corresponding to the quark-diquark structure of a baryon 
is given by 

^ e// = E^ + f + f^ + ^4 2 + ^le|-e(ER)- 7 (E0, (36) 

The Hamiltonian for baryons Eq. (36) looks like the one for mesons [2] 
because we consider basically the string between quark and diquark. There- 
fore the calculations of mean coordinates and electromagnetic polarizabilities 

11 



is the same. But in the case of baryons there are more parameters and the 
analysis is more complicated. 

With the help of Eq. (35) the effective Hamiltonian (4.36) is rewritten as 



n 



e// 



E 



2/^ 



M 



2M 



n + -(RxH) 



7(£xH) 



+ 



2/x 



7T + ^ (^ x H) + 7 (R x H) 



(76|e|-e(ER)-7(E0, (37) 



where (£ x H) fc = e mn k(,mH n . The mass of a baryon J\4(/j,j) is defined here 
as a solution to equation 



Heff* = M^. 



(38) 



In according to the Noether theorem, the momentum II is conserved and 
we can put R = II = in Eq. (37). To find the solution to Eq. (38) we 
can use the substitution 7rfc = —id/d^. Then the mass of a baryon Ai(fij) 
is given by 



M 



where e(/ij, E, H) is the eigenvalue of the equation 



(39) 



<9£ 2 v ^ ' 



7 



+ 2^(exH)' + <7&|e|-7(Efl>* 



e(/i J? E,H)$. 



(40) 



The term (7 2 /(2M)) (^ x H) in Eq. (40) is due to the recoil of the string. 
In nonrelativistic models the effect of the recoil was studied in [13] (see also 
[14]). As we neglected the spin of baryons here, there is no interaction of 
spin with the external magnetic field. It is not difficult to take into account 
such interaction. Eq. (40) is like the equation for mesons [2] and therefore 
we write out the solution for the ground state. If E = H = 0, the solution 
to Eq. (40) when the orbital quantum number / = is given by 



iV 
<$>(p,e,<P) = -Ai{p-a{n))Y lm {e^ 

P 



(41) 
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where p = J p\ + p\ + p\, pk = (2/i6<r) 1 / 3 ^, Ai(p — a(n)) is the Airy function, 
Yi m (8, 0) is the spherical function, N is the normalization constant and a(n) 
are the Airy function zeros so that o(l) = 2.3381, o(2) = 4.0879, o(3) = 5.52 
and so on [15] . For the ground state of baryons the principal quantum number 
n — 1. For the excited states it is necessary to choose the corresponding value 
of n — n r + 1 + 1, where n r is the radial quantum number. The eigenvalue of 
Eq. (40) (at E = H = 0) is 



-1/3 



e(pj) = (2/jp /d (bay l6 a(n) = (<r) 2/i, o(n) 



2/3 



M 



2/i3 (a*i + A^) 



1/3 



(42) 



The condition of the minimum of the baryon mass (39) (dAi(pj) / pj = 0) at 
mi = ?7i2 (and /xi = /J2), with the help of Eq. (41) gives the dynamical mass 
of a diquark: 



A-f 



„<»> + ^> 



a(n) 



3/4 



(43) 



This value is different from one [3] obtained for large angular momentum. 
Using Eq. (43) we arrive from Eq. (39) at the expression for the mass of a 
baryon (see also [3]): 



M(pj 



m 



2/4 0) 



4m? +4M° ) . 



(44) 



To estimate the baryon mass, the value of the string tension a = 0.15 
GeV 2 will be used [3,4]. Neglecting the small current masses of quarks rrij 



(0) 



320 



we find from Eqs. (43), (44) the mass of a diquark for n — 1: p 3 
MeV and the nucleon mass : Ai(pj) = 1.28 GeV [5]. This value of the 
nucleon mass is a little greater then real nucleon mass because spin-spin and 
spin-orbit forces were omitted. 

Now we consider the mean relative coordinates of nucleons on the basis 
of the virial theorem which gives (as for the case of mesons [2]) the mean 
potential energy (U) = 2(T), where (T) is the mean kinetic energy. Then 
using the relation (T) + (U) = e(pj), we arrive at 



(U) = ^e(pj) = |(2 A i)- 1 /3( 6 ^)2/8 o(n) . 



(45) 



Comparing Eq. (45) with the relation (U) = ba{\/£^} gives the following 
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expression 

<\/0 = l^ba)-^a{n) = a~^ ^-J^ . (46) 

In accordance with Eq. (31) the size of the nucleon is characterized by 
the value | z® — (1/2) (z' 1 ) + tS 2 ') |~| b£ |. Introducing the notation r =&£, 
from Eq. (46) we have 

(47) 

The same expression was found in [1,2] for mesons. For the quark-diquark 
system the string tension coincides with those of mesons and therefore the 
quark-diquark system has approximately the same size as mesons. Using 
a — 0.15 GeV 2 and o(l) = 2.2281 we find the mean size of the nucleons 

(■s/p2) = 0.84 fm. (48) 

The experimental value of the charge radii of the proton and neutron are 




rl) = 0.86 fm [16], (r 2 ) = -0.1 13 ±0.005 fm 2 [17]. In accordance with Eq. 



(4.43) the center of mass of the quark-diquark system is situated in the center 
between quark q^ and diquark (q^\ q^) and the mean radius of a nucleon 

is (1/2) (Vr) = 0.42 fm which is the reasonable value. It should be noted 
that charge radii of hadrons are defined from electromagnetic formfactors. 

Now let us consider the case when H = 0, E ^ 0. It is possible to assume 
as an approximation that E || £ [5], i.e. external electric field is parallel to 
the string which connects quark q^ with diquark (q^, q^)- So we neglect 
the rotation of the string. It is justified only for the ground state when the 
orbital quantum number 1 = 0. Introducing the effective string tension 

a eff =*-1e, (49) 

we arrive from Eq. (40) at the eigenvalue e(^,E) = (2/i)~ ' (ba e ff) a(n). 
From (43), (44), by neglecting the small terms containing the current masses 
we find the mass of a baryon 



_M(/x„E) = 4 v /^77 
14 



a[n) 



3/4 



(50) 



which depends on the external electric field. Inserting Eq. (49) into Eq. (50) 
and expanding it in a small electric field one yields 



M(^,E) 



a[n) 



3 / 4 / 2 \ 

U^-^E- q -a-^ E A, (51) 



V v°~ 



where Eq. (43) was used and q = e\ + e2 — e 3 . We write here only terms of the 
expansion Ai (fij, E) in small electric field up to E 2 . The first term in Eq. (51) 
gives the mass of a baryon. The second one is connected with the potential 
energy of a dipole moment of quarks in the external electric field U = — dE. 
From Eq. (33) when the center of mass coordinate R = and | £ | ^>| r\ |, 
the dipole moment of quark-diquark system is d ~ 7^ = (7/6)1". Comparing 
the potential energy of a dipole U = —( , ~f/b)rE (at E || r) with the second 
term of Eq. (51): — (q/y/a)[a(n)/3] 3 ' 4 'E we arrive at the expression for the 
mean relative coordinate r = (2/ A /o : )[a(n)/3] 3/ ' 4 which coincides with Eq. 
(47). Here we define the electric dipole moment of quark-diquark system 
more precisely as compared with the letter [5] and as a result the mean 
relative coordinate of a baryon Eq. (47) coincides with that for a meson. 
The third term in Eq. (51) describes the potential energy due to the electric 
polarizability of a baryon. The polarization potential is given by (see [18,14]) 

U(a,j3) = -±aE?-±j3H 2 , (52) 

where a, /3 are electric and magnetic static polarizabilities of hadrons, re- 
spectively. From Eq. (51) by comparing the quadratic term in E with Eq. 
(52) we arrive at the electric polarizability of a baryon: 



q 2 



a[n) 



3/4 

a' 3 / 2 . (53) 



Let us consider the estimation of the electric polarizability for the proton 
p = uud. There are two possibilities for a proton as a quark-diquark system: 

a) the quark q^ = d and diquark (q^q^) = (uu), so the electric charges 
ei = t2 = (2e)/3, e 3 = — e/3 and parameter q = ei + e2 — e^ = (5e)/3; 

b) the quark q^ = u, diquark (g^gt 2 )) = (ud) and the electric charges 
t\ = e 3 = (2e)/3, t2 = —e/3 and parameter q = e\ + e2 — e 3 = —e/3. It 
should be noted that there are no permutations of quarks here which occur in 
the Green function Eq. (11). Using the value of the string tension a = 0.15 
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GeV 2 and o(l) = 2.3381 from Eq. (53) we find the static polarizability of a 
proton in Gaussian units for two cases 



a) a 



p 



5.56 x 1(T 4 fm 3 b) a p = 0.22 x 10" 4 fm 3 . (54) 



The generalized electric polarizability which is extracted from measure- 
ments of the Compton scattering cross sections is given by [14] 

a = a + Aa, (55) 

e 2 r| e 2 (K 2 + l) 

Aa = 3M^ + -am!T' (56) 

where M^ is the mass of a nucleon, r^ is electric radius and magnetic moment 
of the hadron \x = (e/c)/(2Mw). Using the experimental values of the electric 
radius and magnetic moment of a proton we have Aa p = (4.5 ±0.1) x 10 -4 
fm 3 [14]. From (54) the total electric polarizability of a proton is becomes 



a) a 



p 



10 x 10~ 4 fm 3 b) a p = 4.7 x 10~ 4 fm 3 . (57) 



The configuration of a proton a), when the quark q^ = d and diquark 
(gWgC 2 )) — ( MM ) j s m ore favorable as the experimental values of electric 
polarizability are 

«; xp = (12.1 ± 0.8 ± 0.5) x 10~ 4 fm 3 [19], 

a c ^ = (10.6 ± 1.2 ± 1.0) x 10~ 4 fm 3 [20], 

a^= (9.8 ±0.4 ±1.1) x 10" 4 fm 3 [21]. 

So in the case a) we have a good agreement with experimental data. 

For the neutron n = udd there are also two possibilities: a) the quark 
g(3) _ u anc j diquark (gWgt 2 )) = (dd), so the electric charges t\ = 62 = — e/3, 
e 3 = (2e)/3 and parameter q = e\ ±e 2 — e 3 = — (4e)/3; b) the quark q^ = d, 
diquark (q^q^) = (ud) and the electric charges 63 = e 2 = —e/3, ei = (2e)/3 
and parameter q = e\ ± e 2 — e 3 = (2e)/3. Inserting these parameters into 
Eq. (53) one gives 

a) a n = 3.56 x 10~ 4 fm 3 b) a n = 0.89 x 10" 4 fm 3 . (58) 

For the neutron Aa n = 0.62 x 10~ 4 fm 3 [14] and the generalized electric 
polarizability of a neutron is given by 

a) a n = 4.2 x 10" 4 fm 3 b) a n = 1.5 x 10~ 4 fm 3 . (59) 
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The experimental situation for a neutron is more complicated as there 
are different experimental data: 

a e * p = (12.6 ± 1.5 ± 2.0) x 1CT 4 fm 3 [22], 

a n = (8.8 ±2.4 ±3) x 1(T 4 fm 3 [23]. 

The recent data [23] are close to case a) in Eqs. (59) with quark g® = u 
and diquark (q^q^) = (dd). The value of a n in the situation a) in Eq. 
(59) is close to one obtained in the oscillator nonrelativistic quark model 
[24-26,18,14]. The results of calculations in the framework of the dispersion 
sum rule [27] and the chiral perturbation theory (CHPT) [28] are closer to 
the case a). 

5 Diamagnetic polarizabilities of nucleons 

In according to Eq. (52) for calculating the magnetic polarizability of nucle- 
ons one needs to compare it with Eq. (37). The effective Hamiltonian (37) 
(at E = 0, R = II = 0) can be cast into 

^e// = T~Lo + 'Hinti 

m 2 M Id 2 , . _ 

H —^ HL+ % + ^) Kxii)]2 - (6o) 

where L k = —iek mn C,md n is the angular momentum and d n =d/d£ n . Without 
loss of generality we can choose the direction of the magnetic field on the 
third axis, i.e. H = (0,0, if). Then the Hamiltonian of an interaction of 
quarks with the magnetic field is given by 

Hm( = _« L3 + ^g + ^)(^). (61) 

where L 3 = % {^d\ — £1^2). Considering the small external magnetic field, 
the perturbative theory can be applied. Using the perturbative method [28] 
one arrives at the shift of the energy 



A£ n = (n 



^3+^ + JiWw 



2/x V 8 ^ 2M , 
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n) 



|(n'|-(A#L 3 )/(2^)|n)| 2 
+ Z. ^TT7 > ( 62 ) 

where d is the angle between coordinate £ and magnetic field H. If the first 
term in Eq. (62) is not equal to zero then the second and third terms are 
smaller, and the main contribution to the energy comes from the first term. 
This occurs when the orbital quantum number I > 0. In the case / = 0, the 
shift of the energy due to the interaction with the magnetic field is defined 
by the second term in Eq. (62). 

After averaging Eq. (62) and taking into account the equation (l/47r) / sin 2 fidQ 
2/3 we find for the ground state when / = 0, L3 | 0) = 0, the shift of energy 

a M!4)t« 2 >' (63) 

where (£ 2 ) = (0 | £ 2 | 0), and | 0) means the wave function of the ground 
s-state. We took into account that the first and the third terms in Eq. (62) 
equal zero because L3 | 0) = 0. Here we ignore the spin interactions of 
a baryon with the external magnetic field and therefore only diamagnetic 
polarizability can be defined from Eq. (63). Using the definition of the 
relative coordinate r = 5£ and comparing Eq. (63) with Eq. (52) gives the 
diamagnetic polarizability of a baryon 

^ = -(J + 9^- (64 » 

As required, the diamagnetic polarizability is negative and Eq. (64) is like 
the Langevin formula for the magnetic susceptibility of atoms. The similar 
expression was derived in [2] for mesons. Take into account Eqs. (26), (43), 
expression (64) is rewritten as 



pdia = _ 1_ + 2 V_i (65) 

1 \ 4 H J 6M ' K ' 

where M — /j,[ +/i 2 +/4 = 2/4 = 2y/a [a(n)/3] . To calculate the value 
of (3 dm for nucleons we can use the theoretical magnitude of the mean-squared 
relative coordinate (r 2 ) or experimental data for the size of a nucleon. The 
first way is preferable. Taking into account Eq. (47), the value of M and 
using the approximate relation ( r 2 ) ~ (vr 2 ) 2 , Eq. (65) transforms into 



R dia _ e 2 + V 
1 12a 3 / 2 



a[n) 



3/4 

(66) 



For a proton with the more favorable configuration a), when the quark 
q^ = d and diquark (q^q^) = (uu), parameter q = (5e)/3 and Eq. (5.66) 
(at a = 0.15 GeV 2 ) takes the value 



pdia = _g x 1Q -4 y m 3_ (g ? ) 

This quantity is greater than that found in the nonrelativistic quark model 
[18,14]. The value of the magnetic polarizability extracted from the low- 
energy Compton experiments [19-22] is 

j3 = (3 para + (3 dia , (68) 

where the paramagnetic polarizability is given by [18,14] 

r = 2 ^ IH^I")l 2 (69) 

Here | 0), ] n) are the ground and excited states of a nucleon, respec- 
tively, m z is the third projection of the magnetic dipole operator. The main 
contribution to the paramagnetic polarizability of a nucleon is due to from 
A(1232) excitation [28] and is given by p ^ ra = (13 ± 3) x 10~ 4 fm 3 . The 
intermediate state A is an isovector excitation and contributes to the proton 
and neutron so that f3^ = fi\. To calculate pp ara in the present approach 
one needs to take into account the interaction of spins of quarks with the 
magnetic field. For estimation of the total magnetic polarizability of a proton 
we use the contribution 0^ ra and in accordance with Eqs. (67), (68) we find 

P P = (5 ± 3) x 10' 4 fm 3 . (70) 

This quantity is close (within two standard deviations) to the experimental 
value [22] which was extracted from the measurements of the cross sections 
of the Compton scattering on hydrogen: /^ xp = (2.9 =F 0.7 =F 0.8) x 10~ 4 fm 3 . 
The value (70) is also in agreement with /^ xp = (2.1 =p 0.8 =F 0.5) x 10~ 4 fm 3 
[19] and with the quantity found in CHPT [28]. The value J3 p is small due to 
the partial cancellation of the positive paramagnetic polarizability /^ ara with 
the negative diamagnetic polarizability (3 p i ' ia . The theoretical evaluation, in- 
terpretation and experimental extraction of the total magnetic polarizability 
J3 P is difficult because it is small. That is why this quantity is model depen- 
dent in different schemes. The approach considered allows us to improve the 
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accuracy by using the perturbation in the spin interaction. The next step is 
to take into account such corrections. 

For a neutron with the configuration a) where the quark q^ = u and 
diquark (q^q^) = (dd), the parameter q = — (4e)/3 and Eq. (65) gives 

(3* ia = -5.4 x 1(T 4 fm 3 . (71) 

Using the paramagnetic polarizability f3 p £ ra = (13 ± 3) x 10~ 4 fm 3 [30] 
for a neutron we get from (68) the generalized magnetic polarizability of a 
neutron 

P n = (7.6 ± 3) x 10" 4 fm 3 (72) 

which is in agreement with the experimental quantities 
(3 n = (6.5 =p 2.4 =f 3) x 1(T 4 fm 3 [23] 

(3™p = (3.2 =p 1.5 =F 2.0) x 10~ 4 fm 3 [22] 

According to Eqs. (57), (59), (70), (72) the sum of nucleon polarizabilities 
in our approach is given by 

a p + (3 p = (15 ± 3) x 10~ 4 fm 3 , 

a n + p n = (11.8 ±3) x 10" 4 fm 3 . (73) 

Values (73) are close to the reliable results found on the basis of the sum 
rule [27]. The interesting feature of the experimental data is that for the 
proton and the neutron the approximate equality a p + J3 P ~ d n + J3 n , and 
inequality a p , a n ^> J3 p , J3 n are valid. The last inequality means that both 
a proton and a neutron behave like electric dipoles. In addition, the large 
positive paramagnetic polarizability of nucleons should be cancelled by the 
diamagnetic contribution. It was shown in the framework of the sum rule 
[27] that in the chiral limit, m n — > 0, both polarizabilities a, f3 diverge as 
0(m~ 1 ), and this was confirmed in CHPT [28]. The same behavior of pion 
polarizabilities was found on the basis of quasiclassical calculations [31] in 
the framework of the instanton vacuum theory. Thus, the leading terms of 
electromagnetic polarizabilities are determined by CSB. As pion degrees of 
freedom play a very important role in the nucleon electric polarizability, one 
needs to describe pions in this scheme. In the approach considered here the 
accuracy of the values of electromagnetic polarizabilities can be improved by 
taking into consideration pion degrees of freedom. 
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6 Conclusion 

The model of baryons as a quark-diquark system is very similar to the ap- 
proach for mesons [2]. The effective action for baryons (21) obtained on the 
basis of the QCD string theory is applied here for calculating the mean size 
and electromagnetic polarizabilities of a proton and neutron which are in rea- 
sonable agreement with the experimental data. The preferable combination 
for a diquark is (uu) for a proton and (dd) for a neutron. In this case the- 
oretical values for electromagnetic polarizabilities are close to experimental 
data. 

It should be noted that the quark-diquark ansatz used here generates an 
electric dipole moment, and a result, there is a contribution to the energy 
linear in the electric field which is like the linear Stark effect. 

As a first step, we made some approximations and model assumptions. 
So, spin interactions of quarks treated as a perturbation were neglected here, 
but we took them into account by using the paramagnetic polarizability of 
a nucleon due to the A contribution. As an approximation, it is justified, 
because in Isgur-Karl model of baryons [32] spin-orbit splitting are much 
smaller than expected from one-gluon-exchange matrix elements (spin forces 
were also discussed in [33]). It was also shown [34,35] that the contributions 
from the Coulomb and spin-spin interactions cancel each other. Nevertheless 
all parameters should be taken from the approach and then compared with 
empirical data. Therefore one needs to calculate the paramagnetic polariz- 
ability of a nucleon in our approach. 

Implying small spin-orbit forces in baryons, we come to K* and quark- 
diquark system correspondence (see [2]). It is possible also to consider excited 
states of baryons at nonzero orbital, /, and principal, n, quantum numbers. 
So, the present approach can be applied for studying any baryons (see also 
[34]). 

The electromagnetic polarizabilities of nucleons found are close to the 
values calculated in the framework of the dispersion sum rule [27] and CHPT 
[28]. 

To improve the accuracy of calculations of electromagnetic polarizabilities 
of nucleons, one should take into account pion degrees of freedom because 
the pion cloud contributes substantially to electromagnetic properties. 

The nucleon polarizabilities were also evaluated on the basis of nucleon 
soliton models in [36-42] . 
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